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Abstract 

The formula of expanding the Abel variety theta function restricted 
to Abel subvariety into theta functions of this subvariety is obtained. 
With the help of this formula the solution of differential equations 
with Jacobi theta functions, restricted on a nonprincipally polarized 
Abel subvariety and their translations are rewritten in terms of the 
theta functions of these subvarieties. This is exemplified by the CKP 
equations, the Bogoyavlenskij system, and the Toda g'2^^-chain. 



1 Introduction 

In this article we show how to find solutions of nonlinear equations in terms 
of theta functions of nonprincipally polarized Abelian varieties. 

Solutions of many known integrable equations are expressed in terms 
of theta functions of Riemann surfaces connected with the equations (the 
Korteweg-de Vries equation, the Kadomtsev-Petviashvili equations, geodesies 
on an ellipsoid, Kovalevskaya's top, etc.; see the survey [1]). In some prob- 
lems where the Riemann surface admits a holomorphic involution, the ar- 
guments of theta functions belong to the Prym variety, a subvariety of the 
Jacobi variety of the Riemann surface, and solutions are therefore expressed 
in terms of the theta functions of the Prym variety. The Prym theta formu- 
las are convenient because the theta functions depend on fewer variables and 
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their qualitative analysis is simpler. The Veselov-Novikov equations and the 
Landau-Lifshits equations are among such equations. In the derivation of 
Prym theta functional formulas for their solutions, the fact is essentially used 
that in these cases the Prym variety is principally polarized and therefore 
has, roughly speaking, a unique theta function. 

At the same time, there are well-known equations for which the whole 
dynamics reduces to nonprincipally polarized Prym subvarieties with the po- 
larization type (1, . . . , 1, 2, . . . , 2). For instance, so are the CKP equations 
(the Kadomtsev-Petviashvili hierarchy of type C) [2], the equations of motion 
of a rigid body around a fixed point in the Newton field with an arbitrary 
(homogeneous) quadratic potential which were integrated by O. I. Bogoy- 
avlenskij [3], some generalized Toda chains [4], and the geodesic flows on 
quadrics and 5*0(4) [5,6]. Also, examples are known of integrable systems 
that can be linearized on Abelian varieties with the polarization types (1,3) 
and (1,6) (the Toda g-a^^-chain and the geodesic flow on SO {4) with a special 
metric; see the survey [7]). 

In Theorem 1, we show that for the space of theta functions on a nonprin- 
cipally polarized Abelian variety we can choose a basis that consists of lifts 
of theta functions with characteristics from an isogenous principally polar- 
ized Abelian variety. By Theorem 2, formulas in the Jacobi theta functions 
restricted to Abelian subvarieties and their translations can be rewritten 
in terms of the theta functions of the subvarieties themselves. This is ex- 
emplified by the CKP equations, the Bogoyavlenskij system, and the Toda 

The author is grateful to I. A. Taimanov who posed the problem for useful 
discussions and apt remarks. 

2 Theta Functions of Nonprincipally Polar- 
ized Abelian Varieties 

Let M — C^/A be an Abehan variety with a Hodge form cv. Here A is 
a lattice in C^. There are real coordinates Xi, . . . , X2g in and a form cu 
cohomologous to cu of the shape 
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where 5j are natural numbers and 5j divides Sj^i. In the complex basis 
corresponding to the coordinates SiXi, . . . ,SgXg, the lattice A is given by 
AZ^ + fiZ^, where A is the diagonal integer matrix with diagonal {6i, . . . ,6g) 
and r2 is a symmetric {g x 5()-matrix with Imfl > 0. The tuple {Si, ... , Sg) is 
referred to as the polarization type of M and is an invariant of the cohomology 
class of the form u. Denote by Cm — {^(-^1-^)} the space of theta functions 
of M with the following periodicity properties: 

e{z + \\A) = e{z\A), AeAZ^ (1) 

0{z + OcjIA) = exp(-7riOjj - 2mZj)e{z\A), (2) 

where ej = (0, . . . , 1, . . . , 0) (unity occupies the jth position). The dimension 
of Cm equals Si ■ . . . ■ 6g (see, for instance, [8, 9]). 

An Abelian variety is principally polarized if A is the identity matrix. 
The theta function with characteristic [a, b] of a principaUy polarized Abelian 
variety M = C^/{Z^ + QZ^ } is defined by the absolutely convergent series 

e[a,b]{z\^l) = ^ exp(7ri((n + a),0(n + a)) + 27ri{{n + a), {z + b))), 

where a,b E C^. The function 9[a, b]{z\Q) possesses the following periodicity 
properties: 

d[a,b]{z + m\il) = cxp{2TTi{a,m))d[a,b]{z\n), 

9[a, b]{z + flm\fl) = exp(— 27ri(6, m) — 7ii{m, flm) — 27Ti{m, z))9[a, b]{z\fl), 

where m e Z^. If the characteristic [a, b] is rational then the theta function 
with this characteristic determines a section of the line bundle over M. 

Theorem 1. The theta functions 9[A-^s,0]{z\n), with e e Z^/AZ^, 
constitute a basis for the space Cm- 

Proof. It follows from (1) that the theta function 9{z\A) e Cm has the 
Fourier series expansion 

9{z\A) = ^ a^exp(27ri(m, A"^^)). 

Property (2) implies the following recurrent relation for the coefficients of 
the series: 

am+Asj = exp(7rifijj + 2'Ki{m, /Sr^Q,ej))am- (3) 
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From (3) we infer that the theta function is determined by the coefficients 
{a„t}, where < < 5^, 1 < s < g. Denote by 9^{z\h) the function that is 
given by the absolutely convergent series 



where — exp(7ri(A~^£, QA~^£)) and a^+Am are found from the recurrent 
formulas (3). The functions 6^ constitute a basis for Cm- The recurrent 
system (3) is solvable explicitly: 

Oe+Am — exp(7ri(m, Qm) + 2m{l^~^e, Jim) + 7ri(A~-^£, QA~-^£)). 

Then 

^^(^lA) = ^ exp(7ri((m + A~^£), Vt{m + A"^£)) + 2m{{m + A~^£), z)) 



The theorem is proven. 

The theta functions 9^ — 9[/^~^e,Q\{z\Q) are hfts of the theta functions 
with characteristics frornjhe principally polarized Abelian variety M under 
the isogeny ^ : M M, ^{z) = z, of degree Si ■ . . . ■ Sg. Observe that 
our construction depends on the choice of the principally polarized Abelian 
variety M isogenous to M. If we choose another isogeny then the lifts of the 
corresponding theta functions with characteristics as well determine a basis 
for Cm- 

We now find expansion for the restriction of a theta function of an Abelian 
variety to an Abelian subvariety in the theta functions of the subvariety. The 
curvature form of the bundle associated with the positive divisor is the Hodge 
form; consequently, every positive divisor on an Abelian variety determines 
the polarization on it. 

Denote by M = C^/jAZ^ + hz^} an Abehan variety and by M C M, 
an Abelian subvariety. Suppose that the intersection of the theta divisor 
of M (the set of zeros of the theta function 9{-\Q)) with M determines the 
polarization type on M. Then the restriction of 9{-\(l) to M 

is a theta function of M. Hence, there is an isomorphism (p : C"/{AZ" + 
nZ"} M, where A is the diagonal matrix with the diagonal {Si, ... , 5„) 
and Q is some symmetric matrix with Imfl > 0, such that 9{if{z)\fl) is a theta 




(4) 



9[A-^e,0]{z\n). 
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function of C/fAZ" + QZ"}. Suppose that </?(z) = ^z, where $ is some 

(n X (7)-matrix and 2;''' = {zi, . . . , z^). Since is a theta function 

of C"/{AZ" + nZ"}, wc have the inchision $A C AZ^ and the equahty 
$0 = OP, where P is some integer (n x g()-matrix. Since 

= exp{-7ri{ej,P^QPej) - 27ri{ej, P^<^>z))e{(f{z)\Q), 

P^^ is the identity (n x n)-matrix and consequently fl = P^QP. 
Theorem 2. The following formula is valid: 

e(^{z) - 7|Q) = Yl c.e[A-^£, 0](^ - P^7l^^), 
£eZ"/AZ" 

where 7"'' = (71, ... , 7^) and 

Ce = exjp{7ii{m,ilm) — 27ri(m, 7) 

+27ri(£, A-IPT7) - 7ri(A-^£, f^A-^£)). 

Proof. In view of 

0{ifi{z + A) - 7|Q) = - 7|Q), A e AZ", 

+ ilej) - 7|Q) = 0{if{z) + nPej - 7|Q) 

= exp(-7ri(Pej-, QPej) - 2m{Pej, <^z - -f))e{<^{z) - -f\Q) 

— exp(— 7ri(ej, ilcj) — 27ri(ej, z — P^ ^))9{Lp[z) — 7|0), 

the function d{ip{z) — 7|i7) is a theta function of the Abchan variety M with 
argument z — P^j and hence it expands in the basis theta functions (4). The 
Fourier series expansion of the function 9{ip{z) — ■ylQ) looks like 

exp(7ri(m, Qm) + 27ri(m, $2; — 7)). 

By (4), the coefficient of exp(27ri(A~^£, z — P^7)) in this expansion equals 
exp(7ri(A~^£, QA~-'^£)). Whence we obtain the formula for c^. The theorem 
is proven. 
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3 The Theorem on Expansion of the Prym 
Theta Function 

Here we confine exposition to the case in which the nonprincipally polarized 
Abelian variety M is a Prym variety. 

Suppose that F is a Riemann surface with an involution cr : F — > F 
having 2{n + 1) fixed points Qo, ■ ■ ■ , Q2n+i, n > 0, and tt : F — > Fq = F/cr is 
a projection. The genus of F equals 2g + n, where g is the genus of Fq. The 
surface F possesses a canonical basis of cycles 

Oi, . . . , Op+n; 5i, . . . , dg, bi, . . . , hgj^n-i bi, ■ ■ ■ ,bg 

such that 

(j{aa) + Cla = 0-{ba) + 6q = 0, 1 < « < (?, 

a{aj) + ttj = (j{bj) + 6j = 0, g + ^<i<g + n. 

To the basis of cycles, there corresponds the canonical basis of Abelian dif- 
ferentials 

Ml, ... , UgJ^ni . . . ,Ug 

satisfying the relations 

a*Ua + Ua — 0, a*Uj + Uj — 0, 1 < Oi < g, g + i<j<g + n. 

The cycles ai, . . . ,ag,bi, . . . ,bg project to the canonical basis of cycles on Fq. 
The surface Fq possesses the canonical basis of Abelian differentials wi, . . . , a;^ 
such that 7r*{ujk) — Uk — Uk, 1 < ^ < 9'- Denote by T the matrix of periods 
of ujk, Tjk = f^^i,^)OJk, < j, k < g. Recall that Xr(„^,)Cjfc = Sjk by the 
choice of the basis cycles and the basis Abelian differentials. Denote by 
J ^ c2p+n^|22s+" + nZ'^9+nj ^^c Jacobi variety of F, where O is the matrix 
of periods of the basis differentials of F, and denote by A : F — > J the 
Abelian mapping with basepoint Qo- The differentials 

uJa = Ua + Ua, = Uj, l<a<g, g + l<j<g + n, 

constitute a basis for the Abelian Prym differentials a*u!k — —oJki 1 < ^ < 
g + n. The involution a indices the involution : J ^ J , 

. . . , Zg, Zg^i, . . . , ^g+n) Zi, . . . , 5p) . . . , Zg, • • • ; ^g+n ) Zi, . . . , Zg 
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The Prym variety is the Abehan subvariety 

Pr = {z e J \ a^{z) = -z] C J. 
Denote by (/? : C^"*""/ A — > Pr the isomorphism 

(p{Zi, . . . , Zg+n) = (l/22;i, . . . , IfiZg, l/2z,,...,l/'2Zg+n), 

where A = AZ^~^" + nZ^+", A is the diagonal matrix with diagonal {2, ... ,2 
, 1, . . . , 1) (n units), and 11 is the symmetric matrix with Imll > whose 
entries are 

ba 

Ujk^ J ujk, 9 + l<j<g + n, l<k<g + n. 

The variety Pr is nonprincipally polarized and the dimension of the space of 
the Prym theta functions equals 2^. 

From Theorems 1 and 2 we derive the following 

Theorem 3. The theta functions 9[A~^£, 0]{z\lV), where e — (ei, ■ ■ ■ ,Eg, 
0, . . . , 0) {n zeros), Sj e {0, 1}, constitute a basis for the space of the Prym 
theta functions. The following expansion in the Prym theta functions is valid: 

e{^{z) -^\n) = J2 cMA-'e, o]{z - 7|n), (5) 

where 7 = (71, . . . , 75+^, 71, . . . , %), 7 = (71 + 71, . . . , 7^ + %, 7^+1 . . . , 7g+„), 

— exp(7ri(m£, Jlme) +7ri(£, 7) — 27ri(m£, 7) — 7ri(A~"^£, nA~''"£)), (6) 
mezs 

and rrie stands for the vector (mi, . . . , rUg, 0, . . . , 0, £1 — mi, . . . ,eg — nig) {n 
zeros) . 

In [10, Proposition 5.5], some formula was obtained that connects theta 
functions with the characteristics of the principally polarized Abelian vari- 
eties J, CV{Z3 + TZs}, and C5+"/{Z5+" + nZs+"}. This formula gives 
another way to derivation of (5) and implies that the constants (in 
Theorem 3) equal 9[e, Qi\{S\2T), where 7 = (71 — 71, ■ ■ ■ , % — 7p) and e — 

2(61, . . ■,£g)- 
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4 Applications 



4.1. The CKP hierarchy. This hierarchy is determined by an infinite 
system of Lax equations 

[dn- Bn,dm- B^] = 0, m, n = 1, 3, 5, . . . , 

in the coefficients of the operators 

n-2 
i=0 

that depend on infinitely many variables x = ti, ts, ts, . . .; moreover, the 
equality B* = —Bn must hold. Here B* is the formal adjoint of B^, d — 
d/dx, and 9„ = d/dtn- There is a pseudo differential operator 

V °° 

k=2 

such that Bk = [L^)^ , where {L^)^ is the differential part of and the 
functions Vfe are expressed in terms of V and its derivatives. The first two 
operators of the hierarchy are 

S3 = + Ivd + -dv, B,^d' + ^vd^ + ^dvd^ + wd + ^ + 

2 4 2 4 2 8 

where 

W = ^[^+ ^Vd'V + ^dVd'V + 15{d'vr + ^V{dVr 
+l^y2QV - ^diV - ^dsd'^V - dd^VdV - ^VdsdV + Sd^dV 

.Zoo Z 

The first equation of the hierarchy is 

d^V = -dW - -d^V - 3VdV. 
5 2 

Express the function V (a solution to the CKP hierarchy) in terms of the 
theta functions of the Prym variety. Let F be a Riemann surface with an in- 
volution (7 : r ^ r having 2{n + 1) fixed points Qq, . . . , Q2n+i, n > 0. Fix 
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a nonspecial divisor D ~ Pi + . . . + P2g+n on V and a polynomial R. Take 
a local parameter kr^ in a neighborhood of Qi such that /co" = —k. The 
one-point Baker-Akhiezer function with spectral data {F, Qi, -D, i?(A;)} 
is a function ip{P), P G F, defined (to within proportionality) by the follow- 
ing properties [11]: 

(1) ip{P) is meromorphic on F\Qi and the set of its poles coincides with D; 

(2) the function ip{P) exp{—R{k)) is analytic in a neighborhood of Qi. 
Theorem A [2] If a nonspecial positive divisor D satisG.es the relation 

D + aDr^CT + 2Qi, 

where Cr is the canonical class on T, then the one-point Baker-Akhiezer 
function constructed from the spectral data {F, Qi, D, k~^, xk + t^k^ -{-t^k^ + 
. . .} is an eigenfunction of the operators of the CKP hierarchy; i.e., 



Recall that Cr is the linear equivalence class of the divisor of zeros and 
poles of some meromorphic 1-form. The condition of the theorem means that 
there is a meromorphic 1-form cui with zeros in D + aD and with a second- 
order pole at Qi. 

Let Us be a Prym differential of the second kind with a unique pole of 
order s + 1 (s is an odd number) at Qi and with the zero a- and a-pcriods. 
Denote by Us the vector {Usi, ■ ■ ■ , Usg, Usg+i, ■ ■ ■ , Usg+n, Usi, . . ., Usg), where 
Usj = cOg. The Baker-Akhiezer function ip{x, ts, ts, . . . ; P) looks like [11] 

p p p 

exp (2mx J u!i-\- 27:its J cu^ + 27rit5 J ijJq + . . ^ 
Qo Qo Qo 

9{A{P) - A{D) - Kr + xU^ + ^3^3 + ^5^/5 + . . . |Q) 
^ e{A{P) - A{D) - Kr\n) 

where Kr = — |A(Cr) is the vector of the Riemann constants and 9{.\^) 
is the theta function of the Jacobi variety of F. From the power series 
expansion of ip in a. neighborhood of Qi we obtain the following 

Corollary 1 (to Theorem A) The finite-gap solutions to CKP look like 

V{x, t3, t5, ...) = 2d^ log 9{xUi + tsUs + ^5^/5 + ■ ■ ■ - tI^^), 
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where 7 = A{D) + Kr- 

From Theorem 3 we derive the next 

Theorem 4 The hnite-gap solutions to CKP are expressed in terms of 
the Pryni theta functions by the formula 

V{x, ts, ^5, ■ ■ ■) = 29" log J2 Ce^[A-^£, 0]{xui + tsUs + isC/s + ■ ■ ■ - 7|n), 



where Us = {^Usi, 2Usg, U^g+i, Usg+„), e = (si, . . . , £g, 0, . . . , 0) (n 
zeros), Sj G {0, 1}, and arc found by (6). 

It would be interesting to explain a solution to CKP in terms of secants 
of Abelian varieties. This was done in [12] for soliton equations integrable 
in theta functions of Jacobi varieties and in [13] for equations integrable in 
theta functions of principally polarized Prym varieties (see also [9]). 

4.2. The problem of rotation of a rigid body. Rotation of a rigid 
body S around a fixed point O G S* in a coordinate system Vi, r2, rotating 
together with S in the Newton field with potential ip is described by the 
generalized Euler equations [3]: 

M = Mxa; + — xq; + — X/3 + — X7, (7) 

a — a X uj, $ — P X cu, 7 = 7x0;, 

C/(a,/3,7) = J p{r)^{{r,a), (r,/3), {r,-f)) dridr2dr3, 
s 

where p(r) is the density of the body 5" at r = (ri,r2,r3); a, (3, 7 is the 
orthonormal basis of the fixed coordinate system; and the vectors M and u; 
of kinetic momentum and angular velocity are connected by the relations 



k=l 



with 7jfe the components of the inertia tensor of the body S in the rotating 
coordinate system: 

fik = J p{r) {^ik ^ r| - rirj^j dridr2dr3. 
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On using the isomorphism between the Lie algebra with vector product 
and the algebra of skew-symmetric (3 x 3)-matrices with the commutator 
product, it was demonstrated in [3] that (7) yields the matrix equations 



0, L = BE'^ + ME + u, Q^u- EI, 



where u and B are some symmetric matrices and E is an arbitrary parame- 
ter. For simplicity, we use the same symbols to denote the skew-symmetric 
matrices and the corresponding vectors under this isomorphism. Denote by 
r the smooth completion of the surface that is defined in by the equation 

det{BE^ + ME + u-wl) =0. 

The Riemann surface F is not hyperelliptic and admits the holomorphic in- 
volution cr : F — > F, a{w,E) = {w,—E). The genus of F equals 4, the 
surface Fq = F/cr is elliptic, and the involution a has 6 fixed points. The 
dimension of the Prym variety equals three, and the components cuKt) of 
the angular velocity are equal to 

Aexp(tej^(^^^^^|^), 

where the constants Al and and the vectors zj, zq e J(F), and U e Pr(F) 
are defined in [3]. Prom Theorem 3 we derive the following 

Theorem 5 The components of the angular velocity are expressed in 
terms of the Prym theta functions by the formula 

ujHt) = cxp(tr^) + c,9[A-'e^, 0]{tU + zjlYl) 

' c'Qe(tU + z\U) + c'ie[A-^ei,0](tU + z\U) ' 

where £i = (1,0,0), A = (2, 1, 1), and the constants c* and c'^ are defined 

by (6). 

4.3. The Toda gi2^^-chain. This chain is described by the Hamiltonian 
system with Hamiltonian 

1 ^ 1 
H = - ^Pi+ exp(g2 - qi) + exp(g3 - ga) + exp -{qi + q2 - ^q^) 

1=1 

which is reduced by a change of coordinates to the system 

Y = CX, = {xiyi, X2y2: x^y^), (8) 
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where X^{t) = {xi{t),X2it),Xsit)), Y^{t) = {y^{t),y2{t),y^{t)), and C is the 
Cartan matrix of the Kac-Moody Lie 512^^-algebra which equals 




System (8) admits the Lax representation 



where 
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4 ' ^ 4 ' 4 

It follows from (9) that the operators and dt+B^ commute; consequently, 
the operator takes the seven-dimensional kernel of the operator dt + 
into itself. Hence, the eigenvalues of are independent of time. Therefore, 
the characteristic polynomial Q{n, A) = det(74^ — XE) of ^4^ is independent 
of time: 



Qifi, A) = X{H^{X, F)(l//x+/x) -A6-i/2(^, Y)X^-H:,{X, Y)X''-H,{X, ¥)). 
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The coefficients Hi[X, Y) = Ci (functions of the components A^) are integrals 
of motion. Denote by F the smooth completion of the curve given in by 
the equation 

Ci (1/^ + /i) = + C2 + C3 A^ + C4. (10) 

The spectral curve T admits the two involutions 

T : (//, A) ^ (1///, A), cr : (//, A) ^ (//, -A). 

The involution r is hyperelliptic with fixed points Pj(l, Aj), 1 < i < 6, where 
Aj is a root of the equation A^ + C2X'^ + CsA^ + C4 — 2ci = 0, and -Pj(— 1, Xj), 
7 < J < 12, where A^ is a root of the equation A^ + C2A^ + c^X^ + C4 + 2ci = 0. 
The involution a has 4 fixed points: i?i(/xi,0), i?2(A*2,0), i?3(0, cxd), and 
i?4(oo, 00), where fii and fi2 are roots of the equation Ci(l//i + ;u) = C4, 
and i?4 are points at infinity of the curve given by equation (10). Finite-gap 
solutions to equation (9) are expressed in terms of the theta function of the 
Jacobi variety J of the curve F (see [14]). The arising winding of the torus J 
covers a two-dimensional Abelian subvariety M with the polarization type 
(1,3) [4,7] which is contained in the Prym variety of the involution a. By 
Theorem 2, solutions to (9) can consequently be expressed in terms of the 
theta function of the Abelian variety M. 
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